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An overview of the Lattice technique for studies of the strong interaction is given. 
Recent results from the UKQCD lattice collaboration are presented. These concentrate 
on spectral quantities calculated using full (i.e. unquenched) QCD. A comparison with 
quenched results is made. Novel methods of extracting spectral properties from two- 
point functions are described. 
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1 Introduction 



Lattice Gauge Theory has been applied to the study of the strong interaction in earnest 
for the last 20 years or so. In that time, it has grown from an fledgling, optimistic area 
of research to a well-developed, mature (and optimistic!) discipline. The reason for this 
optimism is that all the approximations involved in the technique are systematically 
improvable. This means that, given enough time on a computer powerful enough, we 
will have predictions for the QCD bound-state spectrum (for example) with arbitrary 
accuracy. Other approaches of studying the strong interaction, while they may have 
their own advantages, cannot make this claim. 

In this talk I will begin by overviewing the lattice method for obtaining hadronic 
spectral quantities via the calculation of n— point functions. I will outline the caveats 
that exist with current lattice simulations as well as detailing some of the successes of 
the approach. The vexed question of performing "full" QCD calculations (i.e. without 
the quenched approximation) is discussed and I will explain why these simulations are 
even more difficult than they at first seemed. In Sec. [3], I will detail some recent results 
from the UKQCD collaboration, focusing on our recent full QCD work. Sec. |] outlines 
a new and promising approach of obtaining physical results from the lattice by using 
the spectral function representation. This is an exciting area of research, and, if it 
reaches its potential, promises to become a standard approach in the future. I then 
summarise the main points raised in this talk in Sec. [|. 



2 Overview of Lattice Gauge Theory 
2.1 The Method 

I attempt here to describe the method normally used in lattice calculations of the 
hadronic spectrum of QCD. There are many excellent reviews of this topic which cover 
the approach in more detail. [1|] 

The conventional lattice approach to QCD spectrum calculations is performed in the 
Euclidean path integral framework. It requires the calculation of n-point correlation 
functions, G n (t), of hadronic interpolating operators, J, in a background sea of glue 
(and sea quarks in the case of full, i.e. "unquenched" calculations). A Monte Carlo 
approach is used to generate these background configurations with the appropriate 
Boltzmann factor e~ s where S is the Euclidean action. There is an obvious analogy 
between this Lagrangian approach to Lattice Gauge Theory, and statistical mechanics: 
clearly the lattice path integral corresponds to the partition function of statistical 
mechanics. 

Note that in lattice simulations there is freedom to choose different quark fields in the 
interpolating operator J compared to those in the Lagrangian. Hence we are able to 
distinguish between "valence" and "sea" quark masses, m val and m sea . 
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Because the QCD Lagrangian contains fermionic fields that appear quadratically, they 
can be integrated out analytically giving the usual determinant factor, det(p+ m). 
Including this factor is a technical headache because it involves a huge increase in 
computational requirements. The usual way around this is to invoke the quenched 
approximation where two things are done: (i) the fermionic determinant is replaced by 
unity; and (ii) the gauge coupling, g, is rescaled so that the physical predictions for a 
particular test quantity (like the rho mass, for example) is in agreement with its exper- 
imental value. The extent to which the quenched approximation reproduces physical 
predictions for other quantities is a measure of its success. It is a remarkable fact that 
for many spectral quantities, the quenched predictions agree with the experimental 
(i.e. the "full" QCD) values to within 10% (see §). 

Despite the success of the quenched approximation, it is obviously essential to perform 
full QCD calculations in order to study the real world. Furthermore, for some quantities 
(such as the deconfinement temperate ||) the quenched approximation is poor, and for 
others (such as the rj' mass) the quenched approximation fails completely. 

In the calculation of the two-point function, G 2 (t), if the exact operator, J exact for the 
hadron in question was used, then G 2 (t) would contain information on that hadron 
and no others. However, since J exact is not known, G 2 (t) receives contributions from 
all hadronic states which have non-zero overlap with J. It is straightforward to show 
that, in this case, G 2 (t) has the following form 

G 2 (t) = 5>e- M *< 

i 

where the sum is over the hadronic states i, and Mj and are the corresponding 
hadronic mass and overlap. Note that since the calculation is performed in Euclidean 
space-time, the excited states are exponentially suppressed with respect to the fun- 
damental state (i.e. the exponentials have real arguments). This means that G 2 (t) 
asymptotes to the two-point function of the ground state hadron as t — > oo. Obviously 
the parameters of the ground state, in particular the mass M = M , can be extracted 
by fitting G 2 (t) to an exponential, e~ Mt , for t sufficiently large. 

2.2 Caveats 

In this sub-section I explain the caveats that one must apply to any calculation using 
the method described above. The main point to make is that the calculations of the 
hadronic properties are performed with input parameter values which are not those of 
the real world. Specifically, these parameters are 

• valence quark mass(es), m val , (typically m val ^ ^m strange ); 

• sea quark mass(es), m sea , (typically m sea ^ m 3tran9e )\ 

• lattice volume, V, (typically V ^ (2 fm) 3 ); 
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• lattice spacing, a, or, through dimensional transmutation, g , (typically a ?Z .05 fm); 

• number of dynamical fermion flavours, Nf. Quenching corresponds to Nf = 0. 
(Typically, N f = or 2.) 

Therefore the mass, M, which was obtained using the procedure above is not the mass 
of the real world hadron, but the mass of the corresponding hadron in a world where 
the quark masses are the same as those input into the lattice calculation, the volume 
is the finite volume used in the simulation etc. So, strictly speaking, M is a function 
of the above 5 input parameters. The final prediction of the real world value, M expt , 
should be obtained by the following extrapolations: 



yjyjyj m val ^few M eV 

vV m sea -> few MeV 

vW V -> oo 

vW a ^0 (1) 

yj Nf — > "2-" i.e. two light flavours for u, d and one heavier for s 

Note that the limit m sea = oo corresponds to Nf = so the m sea and Nf extrapola- 
tions are not independent. While these extrapolations muddy the water significantly, 
many of them are theoretically well-understood and numerically under control. The 
number of -y/'s which appear next to the extrapolations is an indication of how well the 
extrapolation is under control. 

As an example of the quality of the extrapolations, Fig. [I] shows an extrapolation of 
the nucleon and vector meson masses as a function of the lattice spacing, a, taken from 
M. In these plots, data points obtained with both the Wilson and various improved 
lattice actions (designed to have discretisation errors smaller than 0(a)) are shown. 
Fits to the relevant functional forms are included in the figure and the symbols on the 
left of the plot are the continuum extrapolations (which clearly agreement with each 
other). 



2.3 Successes 

The caveats listed in the previous sub-section do not hinder the success of the Lattice 
technique as a means of obtaining accurate predictions from the strong interaction. To 
give an example of the current status of lattice calculations, Fig. shows the hadronic 
spectrum obtained by the CP-PACS collaboration using the quenched approximation. 
There are two features to note. The error bars in the predictions are tiny ( ^ 3%) - 
which is a clear measure of the success of the lattice technique. Secondly, there is a 
small, but statistically significant discrepancy between the lattice predictions and the 
experimental numbers - which is a signal that unquenching is required in order to make 
further progress. 
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Scaling of Vector Mass for 0(a) Improved Quark Actions 




X Wilson Quarks + Wilson Glue (□ CP-PACS) 

o Improved SW Quarks + Wilson Glue 

O Improved SW Quarks + Improved Glue 
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Figure 1: Comparison of scaling violations in the vector meson mass, computed in 
quenched QCD using the standard Wilson action and several different (9(a)-improved 
actions from |4[]. 

While the quenched lattice calculations have clearly matured with precision estimates 
of many quantities of only a few percent, serious calculations involving full QCD have 
only recently begun. Typically the current errors in these calculations are several times 
that of equivalent quenched results. 

3 Recent Dynamical Results from the UKQCD Col- 
laboration 

In this section I review some of the recent results from the UKQCD Collaboration's 
dynamical simulations. |5], || 

3.1 Dependency on Sea Quark Mass 

As was outlined in Sec. [27| all lattice predictions are functions of V, a, Nf etc. In 
this section, the functional dependency on m sea is discussed. Fig. [3] (taken from 0]) 
shows how the lattice spacing, obtained from the hadronic length scaled, r ~ 0.5 fm, 
depends strongly on the sea quark mass, m sea (here expressed in terms of the hopping 
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Figure 2: CP-PACS results for the light hadron spectrum in quenched QCD obtained 
using the Wilson quark action. 

parameter n sea ). These calculations were performed at fixed bare coupling, g . 

This effect has important consequences. Simulations at a fixed bare coupling, go, and 
with several values of m sea , correspond to different physical volumes, and, furthermore, 
different points on the continuum extrapolation a — > 0. Thus finite volume and 0(a) 
systematics are mixed in this case. 

UKQCD have performed two types of calculation. The first was at fixed go for various 
values of m sea in order to calibrate this effect 0. We then performed more sophisticated 
simulations at several values of the parameter pairs (go, m sea ) which were chosen in or- 
der to maintain fixed lattice spacing a ~ 0.11/m (and therefore also fixed volume) J)J. 
This second calculation utilised the "matching" technology of ||. In both these calcu- 
lations, an improved action was used in order to reduce the effect of 0(a) errors. 



3.2 Results 

There is not space to present full details of UKQCD's recent unquenched calculations. 
I discuss the results of only two quantities, and refer the reader to the original papers 
for full details, i, 1 

One of the benchmark quantities of lattice calculations is the static quark potential. 
Fig. [| shows UKQCD's result for this quantity in units of ro from [H]. It can be noted 
that there is little immediate dependency on m sea in this quantity.[| 

1 A closer look however at the data points close to the origin shows a systematic effect which can 
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Figure 3: The inverse lattice spacing plotted against quark mass (expressed in terms 
of 1/ K, sea ) for different lattice sizes taken from The chiral limit is approximately at 
the left margin of the figure. 

Fig. ^| shows the vector meson mass, My, plotted against the pseudoscalar mass 
squared, Mp S from |J. Again, these quantities are expressed in units of r . Quenched 
data which corresponds to the same lattice spacing as the unquenched simulations 
are included as a comparison. The experimental points corresponding to the strange 
mesons are also plotted. 

In summary, the results of ]5], § indicate that the effects of unquenching are small for 
these values of m sea ^ m strange^ ^his motivates the need to move to more physical 
values of m sea « few GeV in the future. 



4 Lattice Spectral Functions 



As outlined in Sec. |2.1| , the conventional method of determining ground state properties 
from lattice simulations is by fitting exponentials to the tails of n-point functions. 
There have been several attempts at developing other strategies for uncovering spectral 
quantities from lattice data. |§ 10|, [11, 112, |13| These all revolve around the spectral 



function (SF) p(s) which can be defined through 



G 2 (t) 



K(t, s)p(s)ds, 



(2) 



where K(t, s) is the kernel function - typically just e~ st for this work. The SF contains 
much richer information on the channel being considered than just the ground state 
parameters. It also has the advantage that theoretical input can be used to guide 



be interpreted as different runnings of the coupling as function of m s 
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Figure 4: The scaled and normalised potential as a function of distance (from []5J). 

its form for large values of s where perturbation theory is valid. In fact, this was 
the approach taken by || [1^] who used (continuum) perturbation theory to derive a 
functional form for p PT (s). This p PT {s) was then used above a certain threshold of 
energy s and a 5 function used for the ground state following the approach of QCD 
Sum Rules. 

A very new and promising technique takes the marriage of lattice data and spectral 
functions one step further. fH]] This approach uses the lattice data itself to determine 
the SF by inverting eq.(0). This is a very numerically technical approach; it is in fact 
an "ill-posed" problem - G2{t) is known only at a small number of discrete values of 
t, whereas the aim is to determine p(s) for a large number of values of s (ideally for a 



continuous range in s). In |L4| a "maximum entropy method" is employed to overcome 
these hurdles. 



5 Conclusion 



I have given a brief overview of the current state of play of lattice gauge theory calcula- 
tions of the hadronic spectrum. The problems of extrapolating the input parameters of 
any lattice simulation to their physical value is emphasised. In particular I have shown 
that interpreting results from naive unquenched calculations can prove difficult due to 
the dependency of the lattice spacing on m sea . A summary has been given of recent 
unquenched results from the UKQCD collaboration for the static quark potential and 
vector meson mass. Finally the new and interesting method of using spectral functions 
in the analysis of lattice data is discussed. 
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Figure 5: Vector mass plotted against the pseudoscalar mass for several sets of different 
sea quark masses from M The experimental points are denoted by * and the quenched 
results are labelled n sea = 0. 
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